Electron-phonon coupling in graphene antidot lattices: 
an indication of polaronic behavior 
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We study graphene antidot lattices - superlattices of perforations (antidots) in a graphene sheet 
- using a model that accounts for the phonon-modulation of the 7r-electron hopping integrals. We 
calculate the phonon spectra of selected antidot lattices using two different semi-empirical methods. 
Based on the adopted model, we quantify the nature of charge carriers in the system by computing 
the quasiparticle weight due to the electron-phonon interaction for an excess electron in the con- 
duction band. We find a very strong phonon-induced renormalization, with the effective electron 
masses exhibiting nonmonotonic dependence on the superlattice period for a given antidot diameter. 
Our study provides an indication of polaronic behavior and points to the necessity of taking into 
account the inelastic degrees of freedom in future studies of transport in graphene antidot lattices. 

PACS numbers: 63.20.kd, 63.22.-m, 71.38.-k, 73.21. Cd 



Recent years have seen a surge of interest in graphene 

- the two-dimensional form of carbon with atoms or- 
dered in a honeycomb lattice. 1 This material shows ex- 
traordinary properties, such as room-temperature ballis- 
tic transport on a submicron scale and the possibility of 
heavy doping without altering significantly the charge- 
carrier mobility. Yet, the usefulness of pure graphene for 
carbon-based electronics^ is limited as the electron trans- 
mission probability across a potential barrier is always 
unity - regardless of the height and width of the barrier 

- a feature akin to Klein tunneling.— Thus the conductiv- 
ity cannot be altered by a gate voltage, the latter being 
a key property of a field-effect transistor. 

Motivated in part by the compelling need to create a 
band gap in graphene, an extensive research effort is cur- 
rently being dedicated to understanding the electronic 
properties of graphene-based superlattices^ A class of 
such structures, made by perforating a graphene sheet - 
graphene antidot lattices - has recently been proposed^ 
These lattices belong to the family of superhoneycomb 
systems^ and can be obtained by patterning graphene 
monolayers using electron-beam lithography, a method 
which allows feature sizes as small as tens of nanome- 
ters. It is worthwhile to stress, however, that - owing to 
recent advances in nanofabrication 7 - - sub-10 nm antidot 
diameters constitute a realistic near- future prospect. 

The electronic structure of triangular antidot lattices 
has been studied theoretically;^ revealing features such 
as the existence of localized midgap states (flat- and 
quasi- flat bands). In addition, transport properties of 
their square-lattice counterparts have been investigated 
experimentally,- showing a transport gap and weak lo- 
calization corrections to the conductance. 

In the present work, we study the influence of phonons 
on the electronic properties of graphene antidot lat- 
tices. We calculate the phonon spectra of selected anti- 
dot lattices using two independent semi-empirical meth- 



ods. We then describe the electron-phonon (henceforth 
e-ph) interaction based on a model that accounts for 
the modulation of hopping integrals by the lattice dis- 
placements (Peierls-type e-ph coupling). 11 Within this 
model, we quantify the effect of phonons by computing 
the conduction-band quasiparticle spectral weight. We 
show that the phonon-induced renormalization is much 
stronger than in graphene itself, with the effective elec- 
tron masses being typically in the range 3.7 — 5 bare band 
masses. 

The triangular graphene antidot lattice {L, R} with 
circular perforations [see Fig. HJa)] has a hexagonal unit 
cell [Fig. [lib)] with side length La and antidot radius Ra, 
where a = 2.46 A is the lattice constant of graphene. If 
we choose a carbon atom (hereafter C atom) on sublattice 
A as the origin, its nearest neighbors are given by the 
vectors Si = (y/3/2, -1/2) a cc , 62 = (0, 1) a cc , and 63 = 
(-V3/2, -1/2) a cc [Fig. fib)], with a cc = 1.42 A being 
the distance between adjacent C atoms. 

Given the large size of unit cells in the antidot lattices 





FIG. 1: (a) A segment of a triangular graphene antidot lattice 
with circular antidots and basis vectors ai and a2. The lattice 
period is ai = a2 1 = La-\/3. (b) unit cell of an antidot 
lattice, with vectors <5i , 82, and 83 specifying positions of the 
nearest neighbors of a carbon atom on sublattice A. 
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FIG. 2: The conduction-band dispersion e c (k) for the {17, 5} 
antidot lattice. The inset shows the L-dependence (9 < L < 
17) of the conduction bandwidth W c for R = 5. 



that we consider - with N at ~ 300—1600 atoms - the cal- 
culation of the electronic structure and the phonon spec- 
tra in the whole Brillouin zone using ab-initio methods 
based on the density functional theory (DFT) is incon- 
ceivable. Instead, we model the band structure using a 
nearest-neighbor tight-binding Hamiltonian^ 



R,m,6 



R+d 



l+ 50R+d m + H.C. 



(1) 



where R designate the unit cells (N of them), d m (m = 
1, . . . ,iV at ) specify the positions of the C atoms within 
a unit cell, 8 stands for the nearest neighbors of the C 
atom at position R + d m , and t ~ 2.8 eV is the nearest- 
neighbor hopping integral. Within our model, the Bloch 
wave functions corresponding to the energy eigenvalues 
e„(k) (n is the band index) are given by V>«k( r ) = 
E m CE k 0mk(r), where mk (r) = N- 1 / 2 £ R e lk ' R ^(r - 
R — d m ) and <p(r — R — d m ) is the 2p z -orbital of a C 
atom at R + d m . To a good approximation, the overlap 
of the 2p z -orbitals on different atoms can be neglected. 

The accuracy of the tight-binding method in the case of 
antidot lattices^ is corroborated by the recently demon- 
strated good agreement with the DFT results for lattices 
with very small unit cells ^ Given that the underlying 
honeycomb lattice is bipartite, the resulting tight-binding 
energy spectrum has particle-hole symmetry^ - a prop- 
erty not retained in the exact band structured Our cal- 
culations show (see Fig.[5]for an illustration) that the an- 
tidot lattices are extremely narrow-band systems: for in- 
stance, in the {L, 5} family (with 9 < L < 17) of lattices 
the conduction bandwidth W c increases from 0.11 eV to 
0.14 eV (see the inset of Fig. in the {L, 7} family 
(12 < L < 17) it increases from 0.020 eV to 0.035 eV. 
The band gap decreases from 0.74 eV to 0.18 eV in the 
{L, 5} and from 0.30 eV to 0.15 eV in the {L, 7} family. 

The phonon spectrum of graphene was studied exten- 
sively, using either ab-initio methods or effective mod- 
els^ In the present work, we calculate the phonon spec- 
tra of graphene antidot lattices using two independent 
methods that have recently been shown to yield very ac- 
curate results for graphene: the fourth-nearest-neighbor 



force-constant (4NNFC) method, in the parametrization 
of Zimmcrmann et ai, and the valence force field (VFF) 
method of Perebeinos and TersoffJ^ 

We study the phonon spectra for the {L, 5} and {L, 7} 
antidot lattices. In each case, we first find the equilib- 
rium lattice configuration by relaxing the atoms until 
forces on them are smaller than 10 -5 eV/A. We then 
construct the force-constant matrix D m p im ipi(R — R') = 
d 2 E tot /du m p(R,)du m ,p,(R'), where M m/3 (R) are the dis- 
placements (J3 — x, y, z) from the equilibrium position for 
an atom at R+d m , and -Etot the total lattice potential en- 
ergy. The normal- mode frequencies uj\(q) and eigenvec- 
tors v A (q) (A is the phonon branch index) are obtained 
from the eigenvalue problem D(q)v A (q) = MuJ(q)v A (q) 
for the dynamical matrix D(q) = ^ R D(R) e~* qR , with 
M being the C-atom mass. 

The salient feature of the obtained phonon spectra is 
that the highest optical-phonon energy at q = is at 
around 195.3 meV - essentially inherited from graphene 
itself and only weakly dependent on L and R - while the 
lowest optical-phonon energy can be as low as 0.69 meV, 
the case of the {17, 5} antidot lattice. The two methods 
used are compared by computing the phonon density-of- 
states -Dph(w) = N -1 ^ x 6[uj — uj\(q)}, which shows 
their good agreement (Fig. [3]), especially in the low- 
energy part of the phonon spectrum. 

Generally speaking, the dominant mechanism of the 
e-ph coupling in the 7r-electron systems is the phonon- 
modulation of the electronic hopping integrals^— - 
Peierls-type coupling^ The latter forms the basis of the 
Su-Schrieffer-Heeger (SSH) model»±2iIS, We thus adopt 
a model comprising an electron term (H c ), the phonon 

term (H = 1) H ph = Eq.A^tqjf^A^A + 1/2), and a 
Peierls-type e-ph coupling term 



Hep= 2 ^ ( & R+d m +S^R+d m + H.c. 



R.m,6,A 
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FIG. 3: The phonon density-of-states for the {17, 5} anti- 
dot lattice, obtained using the 4NNFC (solid line) and VFF 
(dashed line) methods. 
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where 6 = S/\\S\\ is the unit vector in the direction of S, 
1 e^ R (6L q , A + S q ,A) 



UA,R+d„ 



v 



■<(q) (3) 



is the displacement of an atom at position R + d m due 
to the (optical) phonon branch A, and a = 5.27eV/A the 
e-ph coupling constant describing the linear dependence 
of the 7r-electron hopping integrals upon the C-C bond- 
length modulations^ In the last equation, the phonon 
eigenvectors (q) are normalized such that J2 m v m (l) ' 
[v£(q)]* = N at 8 xx , and v^(-q) = [v* (q)]*^l 
In momentum space, the Hamiltonian of Eq. @ reads 



H, 



y Yl 7n„(k,q)a^ k+q a n , k (6 t _ qJv 



+ &qA>> 



k,q, A. 7i 



(4) 

where k creates an electron in a Bloch state VVik 

(eigenstate of H e ) and 7^ n (k, q) stands for the e-ph in- 
teraction vertex function. It can be shown that the latter 
is given by 7*„(k, q) = U„ A „ (k, q) + W* n (k, q), where 



U„i(k,q) = 



v /8Mw A (q) 



m,<5 



*-[v* +S (q)-v^(q)] 



r//^n,k+q\*^n,k 



/ y--^n. k+q \ * x~y 



k 1 



(5) 



is the contribution due to hopping within a single unit cell 
(the indices m + 8 denote neighbors d m + S of site d m and 
the coefficients C^ ,k originate from the aforementioned 
tight-binding band-structure calculation), while 



tU„ A „(k,q) = 



-7==^£ 5-[e^ a <(q)-v*(q)] 



w f _ — ?(k+q)-a/ / ^y7i,k+Q\*/-m.k ika/z^m-k+qx* /on,kl 

x |_e ^o mi j o m -t- e j o mi j 

(6) 

originates from the hopping between adjacent unit cells. 
The prime in the last sum signifies a summation re- 
stricted to the neighbors d TO + S of site d m that sat- 
isfy the condition d m + S = a + d mi for some mi = 
rni(S), with a = ±ai, ±a2, ±(ai — a2). Unlike the 
more conventional Holstein-type e-ph coupling, which 
is completely momentum- independent, the Peierls-type 
coupling depends on both the electron and phonon mo- 
menta. The momentum-dependence of the vertex func- 
tion is more complicated than that of the standard SSH- 
coupling: while the latter is defined on a monoatomic 
lattice, here we study a lattice with a basis {d TO }. It 
is straightforward to check, however, that for the case 
of a monoatomic lattice (jV at — > 1, C T ^ k — !• S nm ) the 
vertex function takes on the standard SSH dependence 
7(k, q) cx [sin(k • a) — sin ((k + q) • a)]^ 

The overlap of the bare-electron Bloch state a„ k |0) and 
the Bloch state |^ nk ) of the coupled e-ph system defines 
the quasiparticle spectral weight Z n (k) = | (^nkla^klO) | 2 , 



a quantity characterizing the renormalization of the elec- 
tron Green's function by the e-ph interaction. Its inverse 
is given byi^Z-^k) = l-d u Re E tl (k,w)| tl)=en ( k ), where 
S n (k, uS) is the self-energy due to the e-ph coupling for 
an electron in the n-th Bloch band. Using the ordinary 
Rayleigh-Schrodinger perturbation theory, one obtains 



l7l(k,q) 



N ' U! 
q,A 



uj\ - e„(k + q) 



i0^ 



(7) 



where the dispersion of the optical phonons is neglected 
in the denominator of the last expression. In particular, 
the Rayleigh-Schrodinger perturbation theory is known 
to describe the properties of coupled e-ph systems at k = 
better than the self-consistent Born approximation.— 
In what follows, we compute the quasiparticle spectral 
weight due to the e-ph interaction for an electron at the 
bottom (k = 0) of the conduction (n — > c) band. From 
Eq. (0, for this special case we obtain 



^(o) = i + ^E 



| 7c A c (k = 0,q)| 2 



q.A 



[e c (0) - E c (q) 



1 2 ' 



(8) 



a quantity that yields the ratio of the effective (in the 
presence of the e-ph interaction) and the bare band elec- 
tron masses: Z ( T 1 (0) = m e ff/ra*. Based on Eq. (|SJ), with 
the aid of Eqs. (0 and (|6]), we evaluate Z~ X {Q) for the 
{L, 5} (9 < L < 17) and {L, 7} (12 < L < 17) families 
of lattices. These demanding numerical calculations are 
performed via parallelization on multiple processors. 

The calculations show a rather strong phonon-induced 
renormalization (see Table Q} compared to graphene, 19 
where Z = 0.93 (at the Dirac points) or larger^ Im- 
portantly, we find a very good agreement between the 
results obtained using the 4NNFC and VFF phonon spec- 
tra, with the effective electron masses m g from 4NNFC 
being slightly larger in all the cases considered. As can 
be inferred from Table HI for fixed L the renormalization 
is larger for the structures with smaller antidot diame- 
ters, which squares with intuition. Another interesting 
feature that we find is a nonmonotonic L-dependence of 
m e ff for given R, with minima for L = 13 in the {L, 5} 
family and L — 15 in the {L, 7} family of antidot lattices. 

Detailed analysis shows that the low-energy phonons 
(below 30 meV) contribute at most 20 percent of the 
overall spectral weight, while among the high-energy ones 
the largest contributions come from two narrow intervals, 
around 173 meV and 194 meV, respectively. These high- 
energy modes typically provide 75 — 80 percent of the 
spectral weight and their salient feature is that they do 
not involve significant atomic displacements in the vicin- 
ity of the antidot edges. 

The obtained strong mass renormalization m e g/ml — 
3.7 — 5 suggests that the charge carriers in the system ac- 
quire polaronic character. Indeed, it is plausible to have 
polaronic charge carriers in a narrow-band system with 
a strong e-ph coupling - a common situation in organic 
semiconductors;^ 5 compared to the latter, graphene an- 
tidot lattices have yet narrower conduction bands and 
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TABLE I: Calculated inverse quasiparticle weights (electron 
mass renormalization) (0) = m e s/mt for various antidot 
lattices {L,R}, based on the phonon spectra obtained using 
the 4NNFC and VFF methods, respectively. 



L 


Z- x {0) 4NNFC 




VFF 


R=5 


R=7 


R=5 


R=7 


9 


5.046 




4.811 




10 


4.836 




4.613 




11 


4.732 




4.509 




12 


4.681 


4.056 


4.452 


4.030 


13 


4.662 


3.827 


4.450 


3.808 


14 


4.668 


3.739 


4.452 


3.725 


15 


4.684 


3.707 


4.459 


3.699 


16 


4.709 


3.758 


4.475 


3.701 


17 


4.733 


3.795 


4.494 


3.756 



lower dimensionality. Given that the system at hand - 
due to its size and complexity - is out of reach of the 
exact-diagonalization methods, we have utilized a per- 
turbative approach. Thus the obtained results are not 
expected to hold quantitatively, but they should still be 
qualitatively valid. These results underscore the rele- 
vance of phonons in antidot lattices and show that trans- 
port in these systems, unlike in graphene, cannot be 
treated as purely ballistic; i.e., it ought to be modelled 
by taking into account the inelastic degrees of freedom. 

It is appropriate to comment on the robustness of our 
results for realistic antidot lattices, which may contain 



hydrogen (H)-terminated edges. The H-atoms couple 
only to the cr-orbitals, while all the bands close to the 
Fermi energy originate from the 7r-electron states. These 
atoms give rise to a small change of hopping integrals 
near an edge due to the ensuing geometrical relaxation, 
leading to a minor band-gap modification.^ Likewise, the 
influence of a handful of H-atoms on the dynamics of 
realistic antidot lattices is also not expected to be dras- 
tic, since - as shown in the present work - the dominant 
phonon modes involve only very small atomic displace- 
ments in the vicinity of edges. Therefore, while related 
issues certainly merit further investigation, our results 
are expected to be largely unaffected in realistic lattices. 

In summary, we have studied the influence of phonons 
on the electronic properties of graphene antidot lattices. 
We have computed the phonon spectra for representative 
antidot lattices and determined the quasiparticle spec- 
tral weight due to the electron-phonon interaction for an 
electron at the bottom of the conduction band. We have 
shown that the phonon-induced renormalization in these 
narrow-band systems is much stronger than in graphene 
itself, providing an indication of the polaronic nature of 
charge carriers. Our study paves the way for future in- 
vestigations of the effect of phonons on the electronic and 
transport properties of graphene-based superlattices. 
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